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Abstract 

We examine the possibihty of distinguishing a supersymmetric gluino from a Kaluza- 
Klein gluon of universal extra dimensions (UED) at the Large Hadron Collider (LHC). 
We focus on the case when all kinematically allowed tree-level decays of this particle 
are 3-body decays into two jets and a massive daughter (typically weak gaugino or 
Kaluza-Klein weak gauge boson). We show that the shapes of the dijet invariant 
mass distributions differ significantly in the two models, as long as the mass of the 
decaying particle mA_ is substantially larger than the mass of the massive daughter ms- 
We present a simple analysis estimating the number of events needed to distinguish 
between the two models under idealized conditions. For example, for m^/m^ = 10, 
we find the required number of events to be of order several thousand, which should be 
available at the LHC within a few years. This conclusion is confirmed by a parton level 
Monte Carlo study which includes the effects of experimental cuts and the combinatoric 
background. 



1 Introduction 



Very soon, experiments at the Large Hadron Collider (LHC) will begin direct exploration of 
physics at the TeV scale. Strong theoretical arguments suggest that this physics will include 
new particles and forces not present in the Standard Model (SM). Several theoretically 
motivated extensions of the Standard Model at the TeV scale have been proposed. After 
new physics discovery at the LHC, the main task of the experiments will be to determine 
which of the proposed models, if any, is correct. 

Unfortunately, there exists a broad and well-motivated class of SM extensions for which 
this task would be highly non-trivial. In these models, the new TeV-scale particles carry 
a new conserved quantum number, not carried by the SM states. The lightest of the new 
particles is therefore stable. Furthermore, the stable particle interacts weakly, providing a 
very attractive "weakly interacting massive particle" (WIMP) candidate for dark matter 
with relic abundance naturally in the observed range. Models of this class include the 
minimal supersymmetric standard model (MSSM) and a variety of other supersymmetric 
models with conserved R parity. Little Higgs models with T parity (LHT), and models 
with universal extra dimensions (UED) with Kaluza-Klein (KK) parity. All these models 
have a common signature at a hadron collider: pair-production of new states is followed 
by their prompt decay into visible SM states and the lightest new particle, which escapes 
the detector without interactions leading to a "missing transverse energy" signature. If this 
universal signature is observed at the LHC, how does one determine which of these models 
is realized? 

One crucial difference between the MSSM and models such as LHT or UED is the cor- 
relation between spins of the new particles and their gauge charges. In all these models, all 
(or many of) the new states at the TeV scale can be paired up with the known SM particles, 
with particles in the same pair carrying identical gauge charges. However, while in the LHT 
and UED models the two members of the pair have the same spin, in the MSSM and other 
supersymmetric models their spins differ by 1/2. Thus, measuring the spin of the observed 
new particles provides a way to discriminate among models. 

Experimental determination of the spin of a heavy unstable particle with one or more in- 
visible daughter(s) in hadron collider environment is a difficult task. One possibile approach, 
which recently received considerable attention in the literature [1-7], is to use angular cor- 
relations between the observable particles emitted in subsequent steps of a cascade decay, 
which are sensitive to intermediate particle spins. This strategy is promising, but its suc- 
cess depends on the availability of long cascade decay chains, which may or may not occur 
depending on the details of the new physics spectrum. It is worth thinking about other 
possible strategies for spin determination. 

In this paper, we explore the possibility of using 3-body decays of heavy new particles to 
determine their spin. The most interesting example is the 3-body decay of the MSSM gluino 
into a quark- antiquark pair and a weak gaugino, 

g^q + q + X- (1-1) 
In a large part of the MSSM parameter space, this decay has a large branching ratio: this 
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occurs whenever all squarks are heavier than the gluinoo Under the same condition, gluino 
pair-production dominates SUSY signal at the LHC We will argue that the invariant mass 
distribution of the jets produced in reaction (II. ip contains non-trivial information about 
the gluino spin, and can be used to distinguish this process from, for example, its UED 
counterpart, ^ q + q + B^/W^. 

It is important to note that the jet invariant mass distribution we study depends not 
just on the spin of the decaying particle, but also on the helicity structure of the couplings 
which appear in the decay (11.11) . as well as on the masses of the decaying particle, the 
invisible daughter, and the off-shell particles mediating the decay. If all these parameters were 
measured independently, the jet invariant mass distribution would unambiguously determine 
the spin. However, independent determination of many of the relevant parameters will be 
very difficult or impossible at the LHC. In this situation, proving the spin-1/2 nature of the 
decaying particle requires demonstrating that the experimentally observed curve cannot be 
fitted with any of the curves predicted by models with other spin assignments, independently 
of the values of the unknown parameters. This considerably complicates our task. Still, 
interesting information can be extraced. For example, we will show that, even if complete 
ignorance of the decaying and intermediate particle masses is assumed, the jet invariant mass 
distribution allows one to distinguish between the decay (11. ip in the MSSM and its UED 
counterpart (assuming the couplings specified by each model) at the LHC. 

The paper is organized as follows. After setting up our notation and reviewing the basics 
of three-body kinematics in Section 2, we present a simple toy model showing how dijet 
invariant mass distributions from three-body decays can be used to probe the nature of the 
decaying particle and its couplings in Section 3. Section 4 discusses using this observable 
for MSSM/UED discrimination, and contains the main results of the paper. Section 5 
contains the conclusions. Appendix A contains the polarization analysis of the decay — > 
q + q + B^/W^ in UED, which sheds some light on the main features of the dijet invariant 
mass distribution in this case. Appendix B contains a brief review of the Kullback-Leibler 
distance, a statistical measure used in our analysis. 

2 The Setup and Kinematics 

We are interested in three-body decays of the type 



where A and B are TeV-scale particles. The main focus of this paper will be on the case when 
A is the gluino of the MSSM or the KK gluon of the UED model, and i? is a neutalino or 
chargino of the MSSM or a KK electroweak gauge boson of the UED; however the discussion 
in this section applies more generally. We assume that q and q are massless, and denote their 

^The main competing gluino decay channel in this parameter region is a two-body decay g gx, which 
first arises at one-loop level and generically has a partial width comparable to the tree- level decay . The 
gluino decay patterns in this parameter region have been analyzed in detail in Ref. [8]. 



A^ q + q + B, 
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four-momenta by pi and p2 respectively. To describe the kinematics in Lorentz-invariant 
terms, we introduce the "Mandelstam variables" , 



U = {Pi+PbY = {PA-P2f , 

t = {p2+PBf = {pA-pi)\ (2.2) 

of which only two are independent since 

s + t + u = 'm\ + m%. (2.3) 

The allowed ranges for the Mandelstam variables are determined by energy and momentum 
conservation; in particular, 

< s < s'"'^^ = {niA - msf. (2.4) 

We will assume that ps cannot be reconstructed, either because B is unobservable or is 
unstable with all decays containing unobservable daughters. Moreover, since the parton 
center-of-mass frame is unknown, no information is available about the motion of particle A 
in the lab frame. Due to these limitations, the analysis should use observables that can be 
reconstructed purely by measuring the jet four-momenta, and are independent of the velocity 
of A in the lab frame. The only such observable is s, and the object of interest to us is the 
distribution dV/ds. This is given by 

ds QAn^ m\ J Eb-pb ~ V) 

where 



E 



ttia + m 



B 



B 



2mA 



PB = sjEl-ml, (2.6) 

and M. is the invariant matrix element for the decay (12.11) . with the bar denoting the usual 

summation over the final state spins and other quantum numbers and averaging over the 

polarization and other quantum numbers of A^ The quantity can be expressed in 
terms of the variables (12. 2p : substitutions 

t —>■ mi — u — — h TTin (2.7) 

rriA-y vfiA-y 



^This procedure should take into account the polarization of A if it is produced in a polarized state. In 
the examples of this paper, production is dominated by strong interactions and A will always be produced 
unpolarized. 
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Figure 1: Dijet invariant mass distribution for the toy model 1 (blue/dashed) and model 2 
(red/dot-dashed) compared to phase space (black/solid) for M^^/niA = 1-5 and niB/mA = 



should be made in |A1p before performing the integral in Eq. fl2.5p . Notice that Eq. (12.51) 
is valid in the rest frame of the particle A; however, since s is Lorentz-invariant, its Lorentz 
transformation is a trivial overall rescaling by time dilation, and the shape of the distribution 
is unaffected. The strategy we will pursue is to use this shape to extract information about 
the decay matrix element A4, which is in turn determined by the spins and couplings of the 
particles A and B. 

To separate the effects of non-trivial structure of the decay matrix element lirom those 
due merely to kinematics, it will be useful to compare the dijet invariant mass distributions 
predicted by various theories to the "pure phase space" distribution, obtained by setting the 
matrix element to a constant value. From Eq. (12. 5p . the phase space distribution is given by 



as the endpoint is approached. 

3 Chiral Structure in Three-Body Decays: a Toy Model 



To illustrate how the chiral structure of the couplings involved in the decay (12. ip can be 
determined from the dijet invariant mass distribution, consider a situation when the particles 

■^Since we are concerned with the shapes of the dijet invariant mass distributions in various models and not 
their overall normalizations, all distributions appearing on the plots throughout this paper are normalized 



0.1. 




(2.8) 



This distributioii§ is shown by a solid black line in Fig. [H Notice that the phase space 
distribution has an endpoint at s = s^^^, with the asymptotic behavior given by 




to have the same partial width T = 
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Model 1: same helicities Model 2: opposite helicities 



-® ^ X ® ^ 



S = Smax 2 ® 1 V 2 ' ® 1 X 

B B 

Figure 2: Momenta (long arrows) and helicities (short arrows) in the A rest frame for s = 
and s = s™*^^ in the two toy models of section 3. 

A and B are real scalars. Introduce a massive Dirac fermion \E' of mass M* > rriA, and 
consider the following two models: model 1 defined by 

£i = VaA^PlQ + ysB^PRq + h.c. (3.1) 

and model 2 defined by 

£2 = VaA^PlQ + ysB^PLq + h.c. (3.2) 
The matrix element for the decay (12. ip in model 1 is given by 



spin 



while in model 2 it is given by 

J2 \M2\' = 2ylyl {{m\ + ml)tu - m\ml) (j^, + ^7^!^) " (^.4) 

spin ^ * * ^ 

The dijet invariant mass distributions in the two models are shown by the blue/dashed line 
(model 1) and red/dot-dashed line (model 2) in Fig. [H Their strikingly different shapes are 
due to the angular momentum conservation law and to the different helicity structure of 
the couplings. To understand this, consider this decay in the A rest frame. In this frame, 
s = 2£'i£'2(l— cos 612)- When s = 0, the quark and the antiquark travel in the same direction, 
as illustrated in Fig. [2l Since A and B have zero spin, the sum of the quark and antiquark 
helicities must vanish for this kinematics. In model 1, the quark and the antiquark have the 
same helicity, and the decay is forbidden for s = 0; in model 2, it is allowed. In contrast, 
when s = s™*^^, the particle B is at rest, and the quark and the antiquark travel in the 
opposite directions. By angular momentum conservation, their helicities must be equal. In 
model 1, this is the case, and the distribution approaches that of pure phase space in the 
limit s — > s™^^. In model 2, this kinematics is forbidden, the matrix element vanishes at the 
endpoint, and the distribution behaves as dT/ds (x {s — s™ax^3/2_ 
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4 Model Discrimination: SUSY Versus UED 

In this Section, we will show that measuring the shape of the dijet invariant mass distribution 
arising from a three-body decay of a heavy colored particle may allow to determine whether 
the decaying particle is the gluino of the MSSM or the KK gluon of the UED model. We will 
begin by comparing the analytic predictions for the shapes of the two distributions at leading 
order. We will then present a parton-level Monte Carlo study which demonstrates that the 
discriminating power of this analysis persists after the main experimental complications (such 
as the combinatioric background, finite energy resolution of the detector, and cuts imposed 
to suppress SM backgrounds) are taken into account. 



4.1 Gluino decay in the MSSM 

We consider the MSSM in the region of the parameter space where all squarks are heavier 
than the gluino, forbidding the two-body decays g qq. In this situation, gluino decays 
through three-body channels. We study the channel 

9{Pa) q{Pi) + q{p2) + Xi{Pb), (4.1) 

where q and q are light (1st and 2nd generation) quarks, and Xi is the lightest neutralino 
which we assume to be the LSP. (Note that many of our results would continue to hold if 
Xi is replaced with a heavier neutralino or a chargino. The only extra complication in these 
cases would be a possible additional contribution to the combinatoric background from the 
subsequent cascade decay of these particles.) The leading-order Feynman diagrams for the 
process (14.1 1) are shown in Fig. [31 the vertices entering these diagrams are well known (see 
for example Ref. [9]). The spin-summed and averaged matrix element-squared has the form 
(up to an overall normalization constant) 



I -A^ MSSM 

spin 



'2 - ^Cl\^F{s, t, u- Ml,) + \Cn\^F{s, t, u; Mr,) , (4.2) 



where 



(t-M2)2 + (n-M2)2 +^(n-M2)(t-M2)- ^^"^^ 

Here m^, m^, M^* and Mr^, are the masses of the gluino, the neutralino, the squarks qi and 
gij, respectively. In order to keep the analysis general, we will not assume any relationships 
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Figure 4: The Feynman diagrams for the KK gluon three-body decay in UED. 

(such as mSUGRA contraints) among these parameters, and will always work in terms of 
weak-scale masses. We also define 

Cl = Tj;Nu-t^{T'j-Q,)Nn, 

Cr = t^Q.Nn, (4.4) 

where = -1-1/2, T| = —l/2,Qu = +2/3, Qd = —l/3,tw = tanO^, and is the neutralino 
mixing matrix^ in the basis {B, W^, H^, H^)- We have neglected the mixing between the left- 
handed and right-handed squarks, which is expected to be small in the MSSM0 Since up and 
down type quarks are experimentally indistinguishable, the dijet invariant mass distribution 
dr/ds should include both the contributions of up-type and down-type squarks. 

4.2 Decay of the gluon KK mode in the UED model 

The counterpart of the decay (14.11) in the universal extra dimensions (UED) model is the 
decay 

9\pa) q{pi) + q{p2) + B\pb), (4.5) 

where and B^ are the first-level Kaluza-Klein (KK) excitations of the gluon and the 
hypercharge gauge boson, respectively. We ignore the mixing between B^ and the KK mode 
of the field, which is small provided that the radius of the extra dimension is small, 
R <C 1/M]Y, and assume that the B^ is the LTP. As in the MSSM case, the decay (14. 5 p is 
expected to have a substantial branching fraction when all KK quarks Q]^ and Q\ are heavier 
than the KK gluon. Note that in the original UED model [11], the KK modes of all SM 
states were predicted to be closely degenerate in mass around M = 1/i?; it was however later 
understood [12] that kinetic terms localized on the boundaries of the extra dimension can 
produce large mass splittings in the KK spectrum. Since such kinetic terms are consistent 
with all symmetries of the theory, we will assume that they are indeed present, and treat 
the masses of the g^, B^, Q]^ and Q\ fields as free parameters. 

"'Wc assume that N is real. It is always possible to redefine the neutralino fields to achieve this. However 
one should keep in mind that the neutralino eigenmasses may be negative with this choice. 

^Large mixing in the stop sector may be present, and is actually preferred by fine-tuning arguments in 
the MSSM (see, e.g., Ref. [10]). However events with top quarks in the final state are characterized by more 
complicated topologies and can be experimentally distinguished fi'om the events with light quarks that we 
are focussing on here. 
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Figure 5: Dijet invariant mass distribution for the UED (blue/dashed) and the MSSM 
(red/sohd) models, compared to pure phase space (black/dotted) for ML^/rriA = MR^./niA = 
1.5 and mB/mA = 0.1. 



The leading-order Feynman diagrams for the decay (14.51) are shown in Fig. HI (We ignored 
the contribution of the diagrams mediated by with i >2, which are suppressed by the 
larger masses of the higher KK modes.) The relevant couplings have the form 



g^Bl [Y{qn) (iYPrQ^ + Y^q^) QI^PlQI + Y{qji) Q^^^Pm + Y^Ql) Qh'PLq] (4.6) 

where Y[qL) = 1/Q,Y{ur) = +2/3 and Y{dR) = —1/3 are the hyperchargesl^ The spin- 
summed and averaged matrix element-squared has the form (up to an overall normalization 
constant) 

J^IA^UEor = YlG{s,t,u-ML.) + Y^G{s,t,u-MR,), (4.7) 

spin 

where M^^ and M^^, are the masses of the left- and right-handed quark KK modes Q\ and 
Qij, and 



^ ' ' ' '~ {t- M2)2 + (n - M2)2 ^ {t- M^){u - M2) ' 



with 



A"''B 



tu 

h2{s,t,u) = 4s(m^ + m|) ^ — ^ [2s{m\ + m^^) + tu — m\m^^) . (4.9) 



niAirt 



A"''B 



4.3 Model Discrimination: a Simplified Analysis 

Armed with the expressions (14. 2 p and (14.71) . it is straightforward to obtain the dijet invariant 
mass distributions for gluino and KK gluon decays and compare them. For example, the 



^The structure of the couplings between the KK gauge bosons and SM (or KK) quarks are unaffected by 
brane-locahzed kinetic terms as long as these terms are flavor-independent. 
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two distributions for a particular choice of parameters, along with the pure phase space 
distribution, are shown in Fig. [51 While not as strikingly different as the two toy models 
of Section [3], the curves predicted by the MSSM and the UED are clearly distinct. (The 
suppression of the UED distribution compared to phase space around s = and s = Smax can 
be easily understood using angular momentum conservation, as explained in Appendix Rl) 
In this section, we will perform a simplified analysis of the discriminating power of these 
distributions, ignoring experimental complications such as cuts, finite energy resolution, 
combinatoric and SM backgrounds, and systematic errors. We will return to include some 
of these complications in the following section. 

The distrubution in each model depends on a number of parameters, including the mass 
of the mother particle m^i, the mass of the invisible daughter m^, and the masses of in- 
termediate particles: {uL,dL,UR,dfi) in the MSSM case and {U\, D\,U]^, D]^ in the UED 
case. We assume that the partners of the up-type quarks of the first two generations and the 
down-type quarks for all three generations are degenerate, and do not include the diagrams 
with intermediate stops (or KK tops) since they produce tops in the final state. Further- 
more, since the Yukawa couplings for the first two generations are small, it is safe to assume 
that m(ui) = m{dL) in the MSSM and m(f/^) = m(D^) in UED. Since an overall rescaling 
of all masses does not affect the shape of the distribution, we need four dimensionless pa- 
rameters to specify the mass spectrum in each model; we use the particle masses in units of 
mA- Experimentally, these four parameters may be very difficult to obtain independently. 
A direct measurement of the masses of squarks/KK quarks may well be impossible, since 
these particles may be too heavy to be produced on-shell. Also, while it is easy to measure 
niA — TUB (one can use the endpoint of the dijet invariant mass distribution or other simple 
observables such as the effective mass [13] or its variations [14]), it is much more difficult to 
measure and individually [15], which would be required in order to obtain m^/m^. 
In this study, we will conservatively assume no prior knowledge of any of these parameters. 
(Of course, if some independent information about them is available, for example the overall 
mass scale is constrained by production cross section considerations, this information can be 
folded into our analysis, increasing its discriminating power.) In addition to the unknown 
masses, the matrix elements in the MSSM depend on the neutralino mixing matrix elements, 
A'^ii and although only the ratio N11/N12 affects the shape of the distribution. Again, 
this parameter is difficult to measure at the LHC, and we will assume that it is unknown; 
fortunately, the effect of varying it is quite small. 

To quantify the discriminating power of the proposed observable, we use the following 
procedure. We assume that the experimental data is described by the MSSM curve with a 
particular set of parameters|l| We then ask, how many events (assuming statistical errors 
only) would be required to rule out the UED as an explanation of this distribution? To 

^Wc checked that the resuhs of our analysis are approximately independent of which model, MSSM or 
UED, is assumed to be the "true" one. 
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Figure 6: Number of events required to distinguish the MSSM and the UED models based 
on the invariant mass distributions of jets from three-body g/G"^ decays. 

answer this question, we scan over 50000 points in the UED parameter space: 

niB/mA = (0...0.5), M(g[)/mA = (1.05 ... 3.0), 
M{D]i)/mA = (1.05... 3.0), M(f/^)/m^ = (1.05 ... 3.0). (4.10) 

For each point in the scan, we compute the KuUback-Leibler (KL) distance (see Appendix [B|) 
between the UED distribution with the parameters at that point, and the "experimental" 
distribution. We then find the "best-fit UED" point, which is the point that gives the 
smallest KL distance among the scanned sample. Finally, we compute the number of events 
required to rule out the best-fit UED point at a desired confidence level. 

The results of this analysis are shown in Fig. O The MSSM parameters used to generate 
the "data" are: = O.lrriA, rn{uR) = m{dfi) = mji,m{uL) = m{dL) = niL, Nn/Ni-^ = 1. 
The parameters mi and m/j were then scanned between 1.05mA and 2mA, and for each 
point in the scan the procedure described in the previous paragraph was performed. Fig. [6] 
shows the number of events required to rule out the UED interpretation of the signal at the 
99.9% c.l. (In the language of Appendix [Bl this corresponds to R = 1000.) In a typical 
point in the model parameter space, about 6000 events are required. For comparison, the 
pair-production cross section for a 1 TeV gluino at the LHC is about 600 fb, corresponding 
to 12000 gluinos/year at the initial design luminosity of 10 fb~^/year. The number of events 
useful for the measurement studied here depends on the branching ratio of the decay (11. ip . 
Since this branching ratio is generically of order one, we expect 0(10'^) useful events/year at 
the initial stages of the LHC running. Thus, at least under the highly idealized conditions 
of this simplified analysis, this method of model discrimination is quite promising in a wide 
range of reasonable model parameters. 

We checked that the conclusions of this analysis are approximately independent of the 
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Figure 7: Number of events required to distinguish the MSSM and the UED models, as a 
function of ms/mA of the "true" modeL 

value of N11/N12 used to generate the "data". They do, however, depend sensitively on the 
ratio mB/mA. as mB/mj^ grows, the MSSM and UED distributions become more and more 
alike. This is illustrated in Fig. [71 which shows the number of events needed to rule out 
the "wrong" model (assumed to be UED) at the 99.9% c.l., as a function of m^/m^ of the 
"true" model (assumed to be the MSSM with m{uB) = ^{dR} = t^{ul) = ^{dL) = l.Sm^ 
and Nil/ N 12 = !)• The UED scan parameters are the same as in Eq. (14.101) . except that we 
vary m^/m^ = (0 . . . 0.9) in this case. It is clear that the discriminating power of the dijet 
invariant mass distribution falls rapidly (approximately exponentially) with growing ms/'n^^. 
This can be understood as follows. The main feature of the invariant mass distributions that 
allows for model discrimination is the presence of the sharp dip at s = in the UED case. 
According to the Goldstone boson equivalence theorem, if the daughter particle B in the UED 
case is highly boosted, the decays into its longitudinal component will dominate. The particle 
B is highly boosted in the vicintiy of s = 0, provided that the mass ratio mB/mA is small; 
as mB/mji grows, the boost becomes less pronounced and the decays into the longitudinal 
component of B are less dominant. This is illustrated in Fig. [8], which compares the ratio of 
partial decay rates into the longitudinal and transverse modes of B for mB/fnA = 0.1 and 
i^b/i^a = 0.5. However, it is exactly the decays into the longitudinal mode of B that are 
mainly responsible for the characteristic dip at s = 0; this feature is far less pronounced for 
the decays into transverse modes. This means that as mg/m^ is increased, the dip gradually 
disappears, and the discriminating power of our observable fades away. 

4.4 Model Discrimination: a Test-Case Monte Carlo Study 

Given the large number of simplifying assumptions made in the analysis of the previous 
section, a skeptical reader may well wonder how meaningful the results presented above are. 
In this section, we will repeat the analysis in a more realistic setting: effects of experimental 
cuts and combinatoric background will be included. We will also bin the distributions, to 
approximate the effects of finite jet energy resolution. Since this analysis involves generating 
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Figure 8: Ratio of the decay distributions of A into the longitudinal component of B to 
the decay distributions into the transverse components of B for mB/raA = 0.1 (solid) and 
'TfiBlrriA = 0.5(dashed). For low rriB/mA the daughter particle is highly boosted at s = 
and will mainly be longitudinally polarized. As increases, the transeverse polarization 
becomes more important. 

large samples of Monte Carlo (MC) events for each model, we were not able to perform a scan 
over the model parameter space, as we did in the previous section. Instead, we will present 
a test case, comparing the MSSM distribution for a single point in the MSSM parameter 
space with the distribution generated by the "best-fit" UED model for that point. 

The chosen MSSM point has the following parameters: = 1 TeV, = OArriA = 100 
GeV, M{Ql) = M{ur) = M{dn) = 1.5 TeV. The corresponding "best-fit" UED point, found 
by the procedure described in the previous section, has the following parameters: rriA = 1.06 
TeV, TUB = 0.15mA = 160 GeV, M{Ql) = M(m)j) = M(ci)j) = 1.6 TeV. (Note that the value 
of mA — TTiB, which can be determined independently, is the same for these two points.) Using 
MadGraph/MadEvent v4 . 1 [16] event generator, we have simulated a statistically significant 
sample (about 20000) of parton-level Monte Carlo events for each model in pp collisions at 
a/s = 14 TeV. The simulated processes are 

pp^qqqqXiXi (4-11) 

in the MSSM, and its counterpart, 

pp qqqqB^B^, (4.12) 

in UED. With the chosen model parameters, the dominant contribution to the processes (14.111) 
and (14.121) comes from pair-production of g/G^, followed by the three-body decay (II. ip . 
which is of primary interest to us. In the MC simulation, we did not demand that the g/G^ 
be on-shell; the full tree-level matrix elements for the 2 — > 6 reactions (14.111) and (14.121) 
were simulated, so that the subdominant contributions with off-shell g/G^ are included. We 
imposed the following set of cuts on the generated events: 

rii < 4.0; AR{i,j) > 0.4; pT,i > 100 GeV; > 100 GeV, (4.13) 
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Figure 9: Left panel: Dijet invariant mass distributions from the MSSM reaction pp — >■ 
qqqqXiXi (green/light-gray), and its UED counterpart pp — >• qqqqB^B^ (blue/dark-gray), 
including realistic experimental cuts and the combinatoric background (Monte Carlo simu- 
lation). Right panel: Theoretical dijet invariant mass distributions ifiom a single gluino/KK 
gluon decay with the same model parameters and no experimental cuts. 



where z = 1...4, j = z + 1...4 label the four (anti) quarks in each event. The first three cuts 
are standard for all LHC analyses, reflecting the finite detector coverage, separation required 
to define jets, and the need to suppress the large QCD background of soft jets. The ^ cut 
is common to all searches for models where new physics events are characterized by large 
missing transverse energy, such as the MSSM and UED models under consideration. De- 
tailed studies have shown that this cut is quite effective in suppressing the SM backgrounds, 
including both the physical background, 4j + Z, Z —>■ uu, and a variety of instrumental 
backgrounds (see, for example, the CMS study [17]). While we have not performed an in- 
dependent analysis of the SM backgrounds, based on previous work we expect that, with a 
sufficiently restrictive ^ cut, one will be able to obtain a large sample of new physics events 
with no significant SM contamination. 

The dijet invariant mass distruibutions obtained from the MSSM and UED MC samples 
are shown in Fig. [91 The distributions are normalized to have the same total number of 
events, since the overall normalization is subject to large systematic uncertainties and we 
do not use any normalization information in our study. Note that for each MC event, we 
include all 6 possible jet pairings; 4 out of these correspond to combining jets that do not 
come from the same decay, and thus do not follow the theoretical distributions computed 
above. In Fig. [9|, we selectecjfl the jet pairs with s < {uia — rriBY; all pairs with larger 

^This selection can be implemented in a realistic experimental situation because niA — mB can be measured 
independently. 
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values of s arise /^from the wrong jet pairings. However, some of the wrong jet pairs do 
have s in the selected range, forming a combinatoric background to the distribution we 
want to study. Nevertheless, it is clear from Fig. [9] that even after realistic cuts ( I4.13P and 
the combinatoric background are included, the distributions in the two models retain their 
essential shape difference expected from the simplified theoretical analysis of the previous 
section. Assuming that the experimental data is described by the MSSM histogram and 
ignoring systematic uncertainties, we find (using the standard test) that about 750 events 
would be required to rule out the UED curve at the 99.9% c.l. Note that this number 
is smaller than those obtained in the previous section, indicating that the performed cuts 
actually enhance the difference between the MSSM and UED distributions. On the other 
hand, the actual discriminating power of the analysis is likely to be somewhat lower than this 
estimate, since the systematic uncertainty in the cut efficiencies was not taken into account 
here. 

Our parton-level analysis does not explicitly take into account the smearing effect due 
to the finite jet energy and direction resolution of a real detector. The hadronic calorimeter 
energy resolution for a jet of energy E can be approximated by 

^.0,05. (l^f. (4,14) 

and is in the 5 — 15% range for the jets that pass the cuts fl4.13p . We can crudely estimate 6s/ s 
to be of order 26E/E, evaluated at E = ^/s. The fractional uncertainty of the measurement 
of s in our analysis is then roughly between 10% (for points with s ~ ^max-j (for 
points with low s). The bin size used in Fig. IHlis of the order of this uncertainty for large s, 
and larger for small s, so we expect that the smearing introduced by binning in our analysis 
provides a reasonable, if crude, description of the expected smearing due to finite jet energy 
resolution. A more detailed investigation of this effect, and other potential detector effects, 
would be required to fully understand the applicability of the proposed method in a realistic 
experimental situation. 



5 Conclusions 

In this paper, we have investigated how the dijet invariant mass distributions ihom. three- 
body decays of a color-octed TeV-scale new particle, such as the gluino of the MSSM and 
the KK gluon of the UED model, can be used to determine the nature of this particle. The 
production cross section for the color-octet state at the LHC is expected to be large, and the 
branching ratio for the three-body decays is significant whenever all squarks/KK quarks are 
heavier than the gluino/KK gluon. If this is the case, the dijet invariant mass distribution 
can be determined accurately at the LHC. The main complication of the analysis is that 
the distributions in the two models we considered depend on a number of parameters in 
addition to the spin of the decaying particle. However, even allowing for complete ignorance 
of these parameters, we found the dijet invariant mass to be a very promising tool for model 
discrimination. 
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The simplified analysis of this paper did not take into account a number of potentially 
important effects. Since the particles involved are colored, the QCD loop corrections to the 
decay amplitudes are expected to be significant, and may modify the tree- level distributions 
we studied. Also, our analysis is performed at the parton level and does not include detector 
effects. While we expect that many systematic effects would cancel out since the analysis 
relies only on the shapes of the distributions and is insensitive to the overall normaliza- 
tion, a better understanding of the systematics is required. We believe that the promising 
conclusions of this preliminary analysis motivate a more detailed study of these issues. 
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A Polarization Analysis of the UED case 

The main feature of the invariant mass distribution of the UED case, which makes it dis- 
tingishable from SUSY, is the dip at s = 0. This feature can be understood by analyzing 
the decay amplitudes of the individual polarization components of the mother and daughter 
particles and considering conservation of angular momentum. As shown with the two toy 
models in Section [HI conservation of angular momentum can lead to suppression of the in- 
variant mass distributions with respect to the pure phase space distribution (12. 8p at s = 0, 
as well as at s = s™^^. The couplings in the UED case have the same chiral structure as the 
second toy model of Section [3l with the quark and antiquark having opposite helicities. The 
added complication in the UED case is that the mother and daughter particles are massive 
spin one particles. We use mz{A) and mz{B) to denote the projections of the A and B 
spins on the direction of the momentum pi of the quark q. These operators have eigenvalues 
■mz{A),mz{B) = — 1,0,+1; the corresponding eigenstates have polarization vectors e_, ei, 
and e+. The transitions among these eigestates are described by a 3 x 3 matrix of decay am- 
plitudes. Using the UED lagrangian (14.61) . we have evaluated these amplitudes and obtained 
the dijet invariant mass distribution corresponding to each entryj^ These distributions, di- 
vided out by the pure phase space distribution ( 12. Sp . are plotted in Fig. [101 At s = the spin 
projections of the quark-antiquark pair sums up to zero, and the final state has no angular 
momentum (see the right panel of Fig. [2]). Therefore the polarizations of A and B must be 
the same. This will result in a suppression of all non-diagonal components in the transition 
matrix at s = 0, resulting in a dip there. At s = s™'^^, however, the spin projections of the 

^For clarity, we only included the contribution of the diagrams with Q\ in the intermediate state. The 
diagrams with Q/j lead to distributions that are identical, up to a parity reflection, to the ones presented 
here. 
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Figure 10: The invariant mass distributions for the decay of individual polarizations, divided 
by the phase space distribution, for m^/m^ = 0.1 and M/mA = 1-5 in arbitrary units. The 
polarization vectors are along the momentum pi of the outgoing quark q. Notice that at 
s = only the diagonal elements are unsuppressed due to angular momentum conservation, 
resulting in a dip of the distribution. 

quark- antiquark pair add up to = +1 (see Fig. [2]). Thus the only allowed decays at Smax 
are the longitudinal component of A to rriz^B) = — 1 and rriz^A) = +1 to the longitudinal 
component of B. Both features at the ends of the distribution can be nicely observed in 
Fig.[IOl 

B The Kullback-Leibler distance 

A convenient measure to quantify how much two continuous probability distributions differ 
from each other is the Kullback-Leibler distance. (For a recent application in the collider 
phenomenology context, see Ref. [4].) In this appendix, we will briefly review this measure. 

Suppose that the data sample consists of events distributed according to the theoretical 
prediction of model T. Consider a second model, S, which predicts a distribution different 
from T. We can quantify the discriminating power of our data sample by the ratio of 
conditional probabilities for S and T to be true, given the data: 

p{S is true|A^ events from T) 
p{T is true I events from T) 
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This equation can be rewritten using Bayes' theorem: 

p{S\N events from T) 



K 



p{T\N events from T) 
p{S)p{N events from T\S) 
p{T)p{N events from T\T) 

where ^(5") and p{T) are the priors - probabihties for S and T to be true before the ex- 
periment at hand is conducted. (In this paper, we assumed that the MSSM and UED are 
a priori equally hkely, so we set p{S) = p(T) = 1.) Suppose that each event i {i = 1 . . . N) 
is characterized by a single variable Sj (in our case, the dijet invariant mass). Since the N 
events are independent, we have 



_ p(s) n..iP(»r'i5) 



For large we can approximate YIn ~ / normalization condition 

Np{s\T) to obtain 



PiS) 



(B.4) 



^^^^exp(-iVKL(T,5)), 
where the KuUback-Leibler distance (also called relative entropy) is defined as 

KL(T,5) := Jd'^^^(j^)p('\'^)- 

It follows that the number of events needed to constrain the probability of model S being 
true, relative to the probability of T being true, to be less than 1/R, is given by 

logfl + logEffl 

KL(r.s)" ' (s-^) 

This number provides a convenient and physically meaningful measure of how different the 
5* and T distributions are. 

Two properties of the Kullback-Leibler distance are worth mentioning in our context. 
First, while this is not manifest from its definition, the KL distance is non-negative, and 
zero if and only if T = S*. Second, it is invariant under transformations s f{s): for 
example, it does not matter whether we consider the jet invariant mass distribution in terms 
of s or rujj = y/s. 
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